Scissors Mode and dichroism in an anisotropic crystal 



in 
o 
o 

(N 

c 



Keisuke Hatada 1 , Kuniko Hayakawa 1,2 , Fabrizio Palumbo 1 
1 INFN Laboratori Nazionali di Frascati, 
c.p. 13, 1-00044 Frascati, Italy 
2 Dipartimento di Biologia, Universita "Padova", 
Via Ugo Bassi 58b, 35121 Padova, Italy 
e-mail: palumbof Olnf . inf n . it 

We suggest that in an anisotropic crystal there should be a new mechanism of dichroism related 
to a scissors mode, a kind of excitation observed in several other many-body systems. Such an effect 
should be found in crystals, amorphous systems and also metallo-proteins. Its signature is a strong 
magnetic dipole transition amplitude, which is a function of the angle between the momentum of 
the photon and the anisotropy axis of the cell. 
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The so-called scissors mode is a collective excitation, 
where two systems can rotate with respect to each other 
conserving their shape. It was predicted to occur in de- 
formed atomic nuclei, where protons and neutrons are 
assumed to form two distinct bodies to be identified with 
the blades of the scissors. Their relative rotation gener- 
ates a magnetic dipole moment which provides the sig- 
nature of the mode. 

In this paper we suggest that a similar excitation 
should exist as well in all systems where there is a de- 
formed atom in a distorted cell, like crystals, for instance 
perovskites, and with the qualifications discussed at the 
end, also amorphous systems and metallo-proteins. In all 
these systems the scissors mode should give rise to a new 
mechanism of dichroism, since the atom can rotate only 
around the anisotropy axis, and therefore cannot absorb 
a photon with momentum perpendicular to it, as shown 
in the Figure. 

The investigation of this mechanism of dichroism has 
several motivations. The most obvious is a direct deter- 
mination of the distortion parameters. In particular for 
metallo-proteins this could provide a new approach to 
the study of their structure. 

A second motivation concerns theoretical calculations 
of dispersive effects because these will be dominated by 
the scissors mode in a channel with its quantum numbers. 

A last reason of general interest is in the comparison 
with similar excitations predicted or observed in other 
systems. 

Indeed, after Richter and coworkers 2 " discovered the 
scissors mode in 156 Gd, it was systematically investigated 
both experimentally 3 and theoretically^ in atomic nuclei. 
Later on it was predicted to occur in metal clusters 5 , 
quantum dots&, Bose-Einstein£ and Fermi § condensates. 
In all these systems, like in the present case, one of the 
blades of the scissors must be identified with a moving 
part (i.e. the valence electrons in metal clusters and 
quantum dots or the atom cloud in Bose condensates) 
and the other one with a structure at rest (the ions in 
metal clusters or the trap in quantum dots and Bosc- 
Einstein condensates). There is one exception in Bosc- 
Einstein condensates, where the simultaneous condensa- 
tion of two different species of atoms has been obtained 9 , 



thus realizing a situation closer to the nuclear case. Until 
now, as far as we know, among all these systems the scis- 
sors excitation has been discovered only in Bose-Einstein 
condensates 10 , and are under experimental investigation 
in metal clusters^. 

Also in these cases the study of the scissors mode pro- 
vides several pieces of distinctive information. In Nuclear 
Physics it is related to the superconductivity of deformed 
nuclei, in Bose-Einstein condensates it provides a signa- 
ture of superfluidity, in metal clusters it is predicted to 
be responsible for paramagnetism. 

In the present case numerical estimates are easy for 
ionic bindings, and as an example we consider the dis- 
torted crystal LaMnC>3, which has a perovskite structure 
and has recently attracted attention due to colossal mag- 
netoresistance effects. We find that the lowest excited 
state has an energy of about 4 eV with a Ml transition 
amplitude of 1.5 a.u. We think that these values can be 
taken as an order of magnitude also for amorphous sys- 
tems and metallo-proteins which have covalent binding. 

We study the motion of the atom in the cell frame 
of reference. We assume the z-axis parallel to the 
anisotropy axis. We do not need to specify the directions 
of the x-, y-axes because we assume axial symmetry. We 
will also use the intrinsic frame of reference of the atom, 
with the £-axis parallel to its symmetry axis. We call 
9 the angle between the £- and the z-axis, and we de- 
note the components of £ in the cell frame according to 
Cx = sin (9 cosy, Cy — sin 9 simp and £ 2 = cos 9. 

The potential energy of the atom depends only on 9, 
and we take it to be 



V = - C sin 2 I 
2 



(1) 



We identify two of the Euler angles of the intrinsic frame 
with tp and 9, the third angle 7 defining the orientation 
of the atom around the £-axis. We assume the kinetic 
energy of the atom to be that of a rigid rotator. 

In quantum mechanics rotations around the symmetry 
axis are forbidden, and therefore the wave functions must 
satisfy the constraint £ • J = — i = 0, where J is the 
angular momentum the atom. Its components in the cell 
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frame become then those of a point particle so that 
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(2) 



where I\ is the moment of inertia along the one intrinsic 
axis. There remains another constraint. Indeed the two 
orientations of the atomic symmetry axis cannot be dis- 
tinguished from one another, so that the admissible wave 
functions must be symmetric under the operation 



9 



7T + If. 



(3) 



Needless to say, the Hamiltonian is symmetric under this 
operation. The scalar product for such a system remains 
defined by the measure of integration over the Euler an- 
gles, which for functions independent of 7 is 

(Vi|V/)=/ dip f d0Bhieil>t(9,<p)1>/(9,<p). (4) 
Jo Jo 

Exploiting the r-symmetry we will write the matrix ele- 
ments in the form 



2tt 



71-/2 



{ipi \ 0\ip f ) = / dip d9 sine 
Jo Jo 

xi>*(6,tp)[O(<p,0)+ 0(Tr + <p,7r-6)]i> f (9, V ). (5) 

Notice that accordingly the wave functions are normal- 
ized to ^. 

The exact eigenfunctions and eigenvalues are knowni^, 
but in the cases of physical interest the quadratic approx- 
imation is sufficient 



Because of the r-symmetry in the region ir/2 < 9 < tt, 
wehavex n ,| m |(7r-0) = (-l) m Xn,\m\(9)- The eigenvalues 
are 

E n<m = ?uj(2n+\m\ + l), {n = 0,1,2...) (10) 
where 



(11) 



The hamiltonian © is exactly the intrinsic hamiltonian 
which appeared in the nuclear model. 

The first excited states have n = 0, m = ±1, and as 
we will see, are the only states strongly coupled to the 
ground state by electromagnetic radiation. As shown in 
the Figure, they describe the precession of the atom at 
an angle 9 ~ Q Q around the axis of the cell. Due to the 
small value of #0 the atom is essentially polarized. This, 
as explained below, has an important consequence on the 
experimental measurements. 

Now we evaluate the Ml and E2 electromagnetic tran- 
sition amplitudes. We will find that the Ml-amplitude 
is large, being proportional to 1/0q, while the E2 one is 
proportional to 9q and therefore small. 

The intrinsic magnetic dipolc moment of the atom is, 



M (Ml , /i) 



^- J drj(r) ■ rAv|ry 1(i (f) 



The atomic electric current 



J O) = — P( r )p( r ) 



(12) 



(13) 
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< 9 < 7T 



1/2C9 2 
1/2C(tt-9) 2 



< 9 < it/2, 

Tt/2 < 9 < TT. 



(6) 



We recognize that this Hamiltonian, in each of the 9- 
rcgions above, is exactly that of a two-dimensional har- 
monic oscillator, provided we identify 9 with the polar 
radius. As shown by the following estimates, the fall off 
of the wave function is so fast that we can extend without 
any appreciable error the integral over 9 up to infinity. 

In the quadratic approximation, in the region < 9 < 
tt/2 the eigenfunctions are 
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where 



Xn,\m\(Q) — u n,\m\ (jf^J 



H exp f_|l^(M>f£ 
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are Laguerre polynomials 12 and 



■m 



1) : 



(7) 



(8) 



(9) 



depends on the electric charge e, the mass m e , the mo- 
mentum p ( r ) of the electrons and the atomic density 
p ( r ) . The intrinsic atomic magnetic moment M. can be 
rewritten 



M ( Ml , n ) = 



2 m P c V 4 7T 



(14) 



where J M is the atomic angular momentum in the intrin- 
sic frame 



J» = \ J drp(r) (r Ap) M . 



(15) 



By the rotation R = e lJ y e e lJzV we can transform it 
to the cell frame 



RJ U R^ 



e'^l icot 



dip' 
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p = 



r de 



p = ±1 



(16) 



Taking in to account the r-symmetry the Ml transition 
amplitude is 



I M.(M1, jLt) I ^0,0 ) : 
d 

x 5 mil 2{xnM I -qq Ixo.o ) 



e h 

2 m e c V 4 7r 



(17) 
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Performing the change of variables £ = 9 2 j9\ 
d 

{Xn,M I qq \Xo,a) = v n ,\v,\ ^0,0 #o 



at; 

we get its final expression 
(ip n , m \M(M l 5j u)]^o,o> = 



4 0) (0 



e ?i 



(18) 



m e c8o V 4-7T 



<W<W (19) 



As anticipated only the states n = 0, m = ±1 are excited. 
The fact that the amplitude vanishes for p = is due to 
the fact that the atom can rotate only around the z-axis. 

Next we evaluate the E2 transition amplitude. The 
intrinsic quadrupole moment is 

Q 2M = ef drp(r)r 2 Y 2 ^(r). 

We parametrize the charge density according to p(r) = 
s(r(9, ip)—r) po where s is the step function and r(0, (p) — 
r (1 + a 2 oY 2 o(r)), so that 

Q2p. — epo^o a 20 <W (20) 
Transforming to the cell frame we have 

M(E2,n) = £)Q a „<2i/|JZ|2M> =Q 2 o4 2 >K* w ,(21) 



where 



'(2-/^)1 



P 2 M (cos0). 



So the -E2 transition amplitude becomes 

( tpn,m | M(E2, p) | "00,0 ) = <W Q20 



(22) 



(2) 



2 S rnfl Q20 



dOe XnM d^(0) Xo ,o, (23) 



since dfli-K-6) = R m c^(0). 
excitation of the lowest state is 



v (2) 



The amplitude for 



Again we find that only the states m = ±1 are excited. 

The excitation energy and the electromagnetic transi- 
tion amplitudes depend only on the parameters C and 
I. To perform an estimate of these parameters we must 
distinguish the ionic from the covalent binding. As a 
prototype of the first one we consider LaMnC>3. In this 
crystal the Mn 3+ ion is surrounded by six O 2- ions sit- 
ting at the vertices of an octahedron which is elongated 
in the z-direction as shown in the Figure 



The restoring force originated by the quadrupole- 
quadrupole electromagnetic interaction is a function of 
the rotation angle 9 



V(6) 



(25) 



where Squ = qu{9 ) — qa(0). This equation holds if the 
charge of the ligands is external to the atom. Qa are 
the diagonal components of the electric quadrupole mo- 
ment density of the cell (we assume Qij = 0, i / j) 
and qu(6) are the electric quadrupole components of the 
atomic electrons after their coordinates are rotated by the 
angle 9 through the s-axis. The expectation values are 
taken with respect to the ligand and atomic wave func- 
tions respectively. To second order in 9 we have Squ = 0, 
<S<?22 = -( - 933 ) 8 2 , Sq 3 3 = ( 922 - <?33 ) # 2 ■ We assume 
that the atomic charge has an axially symmetric ellip- 
soidal shape with axes c = ro ( 1 + 26), a = Tq{1 — 5), 
while the ligands are point-like with coordinates 

(0,0,P (1 + 2 A)), (0,0, -#(,(1+ 2 A)), 

(P (1-A),0,0),(-P (1- A), 0,0), 

(0, R (l - A), 0), (0, -J2o(l - A), 0). (26) 

Therefore the following relations hold 

(q 33 ) = -2(g22) = -2{q u ) 

' ^33) - -2(4^22) = -2(40n> ( 27 ) 



P 5 



P 5 



P 5 



and 



(q 33 )^^^r 2 6, (±g 3 3)^-36 e Z ; A. (28) 

Zi, Z a are the charges of the ligands and the atom re- 
spectively Assuming for simplicity A — 6 we get 



C=^ Zl Z a e 2 ^6*. 



(29) 



Since the moment of inertia of the atom, under the 
simplifying assumption of constant density in the ellipsis, 
is 



2 m e Z a 2 / , r \ 
h - : r Q (1 + 5), 



(30) 



(24) we finally get 



162 Zi e 2 S 2 
m e Rl 



Rk 



18^2 Z a 6^Tl 



(31) 



According toi&, in LaMnO,3 the distortion parameter 
is A = 0.06, while r = 1.42 A, R = 2.02 A, Z t = 
2, Z a = 22. With the above values we get hui — 
4.0 eV, 9 — 0.33 radiants. The small value of #0 jus- 
tifies the extension we did of the integrals over 9 up to 
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(a) (b) 

FIG. 1: In (a) the atom is in the ground state and the incom- 
ing photon has circular polarisation. In (b) the atom carries 
the angular momentum of the absorbed photon and precesses 
at an angle Go around the cell axis generating a magnetic 
moment proportional to l/6o, with a variation of quadrupole 
moment proportional to 9q. 

infinity. The exact valued of ui differs from the above by 
less than 6%. 

A major difference with respect to the nuclear case 
has already been emphasized: In standard nuclear ex- 
periments the nuclei are unpolarized, while in crystals 
the deformed atoms are essentially polarized. As a conse- 
quence the electromagnetic cross sections have an angular 
dependence which should make their measurement possi- 
ble. Indeed the cross section is maximum when the mo- 
mentum of the photon is parallel to the anisotropy axis 
and vanishes when it is perpendicular to it. Therefore 
a differential dichroic absorption with linear polarization 
of the incoming photon parallel and perpendicular to the 
anisotropy axis should be able to single out the magnetic 



dipole transition due to the collective scissors mode, pro- 
vided that around its energy the linear dichroism orig- 
inating from electric dipole (El) transitions be negligi- 
ble or structureless. This latter feature should obviously 
be derived from a realistic description of the continuum 
particle-hole excitations based on the band structure of 
the crystal. In this connection we note that the damping 
of the collective mode is expected either to be zero, if its 
energy falls in the band gap, or rather small if it falls in 
the particle-hole continuum, due to the smallness of the 
matrix clement connecting the mode to the particle-hole 
states. 

In a triaxial deformation an energy splitting of the col- 
lective mode is expected, as predicted in nuclei^. Such a 
splitting might obviously provide additional pieces of in- 
formation, but in a more realistic microscopic description 
based on the Random Phase Approximation, as done in 
Nuclear Physics, one should also expect a fragmentation 
of the mode that might mask it and make its experimen- 
tal detection more difficult. 

As a last remark we may note that inelastic photon 
scattering is another technique able to detect the scissors 
mode. 

The case of covalent binding cannot be analyzed in the 
same way. First, the decoupling of the intrinsic degrees of 
freedom from the collective (rotational) ones, which is at 
the basis of our semiclassical treatment, is less obvious. 
Second, the estimate of the restoring force and of the mo- 
ment of inertia requires a numerical calculation. Indeed 
we expect that the restoring force will result smaller. Yet 
we deem that the scissors mode will not be precluded by 
the nature of the covalent binding, and that the above 
numerical estimates, should still remain valid as an order 
of magnitude. 
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